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Abstract 

Using chiral symmetry we investigate the leading SU(3) violation in the complete set of quark 
twist-3 light-cone distribution functions of the pion, kaon, and eta, including the two-parton distri- 
butions 4> P M , ^) a M , and the three-parton distribution cffi M . It is shown that terms non-analytic in the 
quark masses do not affect the shape, and only appear in the normalization constants. Predictive 
power is retained including the leading analytic m q operators. With the symmetry violating cor- 
rections we derive useful model-independent relations between (/>^' CT ' 3 , <^ ,cr ' 3 , rf^K+^KO, and <j^o\-- 
We also comment on the calculations of the moments of these distributions using lattice QCD and 
light-cone sum rules. 
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I. INTRODUCTION 



Meson light cone distribution functions (LCDFs) play important roles in high energy 
hadronic exclusive processes [1-5]. The same LCDFs contribute in many processes relevant 
to measuring fundamental parameters of the Standard Model [6], such as B — > niu, r\tv which 
give the Cabibbo-Kobayashi-Maskawa (CKM) quark- mixing matrix element \V uh \, B — > Dn 
used for tagging, and B — > nn, Kir, KK, nr], . . . which are important for measuring CP 
violation. 

With the increasing accuracy in data from the B factories, the flavor dependence in 
LCDFs becomes important to understand the flavor symmetry breaking in processes like 
B — > MM' and B — > MV , where M and V are pseudoscalar and vector mesons, respectively. 
In Ref. [7], chiral perturbation theory (ChPT) was first applied to study the leading SU(3) 
symmetry breaking effects in twist- 2 LCDFs It was shown that terms non- analytic 

in the quark masses do not affect the shape, and only appear in the normalization constants. 
Furthermore, with the symmetry violating corrections useful model-independent relations 
between (f>%, <f>^, 4>k+ k°i and <p^ K _ were derived. 

Recently ChPT has also been applied to the computation of hadronic twist-2 matrix 
elements [8, 9]. Many applications have been worked out, e.g., chiral extrapolations of lattice 
data [10-13], generalized parton distributions [14-16], large N c relations among nucleon and 
A-resonance distributions [9], soft pion productions in deeply virtual Compton scattering [17- 
19], pion-photon transition distributions [20] and exclusive semileptonic B decays [21]. The 
method is also generalized to the multi-nucleon case [19, 22]. 

In this paper, we apply ChPT to higher twist matrix elements. We study the leading 
SU(3) symmetry breaking in the complete set of quark twist-3 LCDFs, including the two- 
parton distributions (p p M , (ff M , and the three-parton distribution 0^. Although those twist- 
3 contributions are parametrically suppressed by inverse powers of large scales, they are 
numerically important in B — > MM' , B — > MV [2] and the meson electromagnetic form 
factor [23]. 

In the following sections, we will first summarize our main results on the leading SU(3) 
symmetry breaking of LCDFs, then study the three twist-3 LCDFs sequentially. Finally, we 
will comment on calculations of those quantities using lattice QCD [24-26] and light cone 
sum rules [27-32]. 
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II. SUMMARY OF RESULTS 



The two-parton LCDFs for pseudoscalar meson M are defined by the matrix element of 
the quark bilinear operator [34] 

(M b {p)%{ V -n) [|n,-|n]A a g a (-|n)|0) = ^ j\x e^'^^f^ <f, p M (x, 

-75 (/1AM ~ §/^P^0mM) }^(1) 
where n is a constant light-like vector, n 2 = and our octet matrices are normalized so that 

y y 

tr[A a A 6 ] = 5 ab . [—n, — —n] denotes the Wilson line connecting the quark bilinear located at 
different space-time points on a light cone. 0^- is a twist-2 LCDF while (ff M and <$> a M are 
twist-3 LCDFs for pseudoscalar meson M. x is the quark momentum fraction and \i is the 
perturbative QCD renormalization scale. For simplicity we work in the isospin limit and the 
MS scheme, and normalize the distributions so that fdx (f> l M (x, fj,) — 1 with i = P, p, a. 

Generically from chiral symmetry the leading order SU(3) violation for <p l M takes the form 
[M = n,K, rj] 

2 2 

(f) l M (x,n) = (f) l (x,n) Tj^h E M l ( x ^) ln ( 1 ^f) 

+F^\x,fi,n x )\ . (2) 



The functions 0q , and F^ 1 are independent of m q , and are only functions of Aq CD , 

/i, and x. F^ 1 also depends on the ChPT dimensional regularization parameter which 
cancels the \n(m 2 N / dependence, so by construction <p l M is independent. Throughout 
the text, 7] denotes the purely octet meson. 

Using ChPT, we found very similar results to the twist-2 case [7] at leading order in 
SU(3) violation [<D(m q )\ : 

1) The twist-2 and twist-3 LCDFs are analytic in m q , meaning that 



E% % (x) = 0, Etf(x)=0, Efl(x) = 0. (3) 

The leading logarithmic corrections can all be absorbed by the normalization constants f l M . 
2) By charge conjugation and isospin symmetry, 

fax) = #(1 - x) , <fi,(x) = tfi,(l - x) . (4) 
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And by isospin symmetry 

<P K +(x) = ^ (x) , <f K -(x) = (x) . (5) 

These two equations are true to all orders in m q . 
3) 

<Pk+(x) ~ ^+(1 ~x) = P K -(1 -x)- fa.(x) = (m s - m) 5<P\x) , (6) 

where dtp 1 (x) is m q independent. 

4) Gell-Mann-Okubo-like relations exist among the octet mesons 



+ 30; (x) = 2[^ K+ {x) + p K _{x)\ . (7) 

5) The three-parton LCDFs also have relations similar to l)-4) [see Eqs. (66-68)]. 

6) Statements l)-5) are still true in quenched and partially quenched simulations, and 
with the leading finite volume and finite lattice spacing corrections. 

7) The light cone sum rule results [32] for twist- 3 Gegenbauer moments are consistent 
with the ChPT prediction [i = p, a] 

4a«> 1 = af + 3af . (8) 

The analogous ChPT relation for twist-2 moments puts a tight constraint on the numerical 
values of a^' P . 



III. TWO-PARTON LCDFS 



The operator product expansion of the non-local quark bilinear operator in Eq.(l) gives 
rise to the twist-2 operator Of' a and twist-3 operators O p ^ a and (Ol' a )' lu : 

Of' a = ^ 7s A a (m-lf) k ^ (9) 
CT = ?7sA a (in (10) 
(Orr = ^7feA° (in ■ ^) k+1 ^, (11) 

where D = D — D and k = 0,1,2,... Here having the vector indices dotted into 
n no . . . n m+i) has automatically projected onto the symmetric and traceless part. The matrix 
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elements of these operators yield 

(M b (p) |Of a |0) = -i/jfrPVp)** 1 , (12) 

<M b (p) |OriO> = -i/j^" • p) fc , (13) 

<M b (p) | (0*T |0> = - k -^r M 5 a \n ■ P ) k (p"n" - p"n") (z% , (14) 
where z = 1 — 2x and the moments are defined as 

(z k Y M = [ 1 dx{l-2x) k <P M (x) . (15) 
Jo 

Following similar procedures as in the twist- 2 case [7], we analyze the twist-3 matrix 
elements. 



A. The (j) p M Analysis 

It is convenient to write 

or = 0K* - GEL, (is) 

where 

O p k % R ^ L X a LR [in-V] k ^ R , (17) 

and similarly for O kRL . ^l,r = [(1 T 75)/2]V ; is the left(right)-handed quark field. The 
distinction between \ a LR and \ a RL is only for bookkeeping purposes. We will set \ a LR = 
\ RL = \ a at the end. 

When a = 3 or 8, 0\f transforms simply under charge conjugation (C), being even when 
k is even, and odd when k is odd. The meson states ir° and rj (i.e., M 3 ' 8 ) are C even. Thus 
from Eq. (13), (z h ) P n0 vanishes when k is odd due to C (and using isospin the same applies 
for M = 7r ± ). For all a's the operator would transform as 

C~ 1 O p k ' a C= (-l) k O p k ' a ; (18) 
if we demanded that, under the C transformation 

To construct the corresponding hadronic ChPT operators, we define E = exp(2i7r a A a //) 
and m q = diag(m, m, m s ) = m q . Under a chiral SU(3)l x SU(3) r transformation, we have 

E -> LHR\ m q -> Lrn^i? 1 " (20) 

A£ fl - L\ a LR R\ \ a RL -> i?A« L Lt. (21) 
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Under charge conjugation, E — > , m q — > mj, while X%r,rl transform according to 
Eq. (19). At next to leading order (NLO) in the P 2 /A^ and m 2 M / k 2 x chiral expansion, 



(22) 



where the C's are the leading order (LO) and the C's are NLO. The sum on i runs over 
hadronic operators which have the same transformation properties as Ol' a . The ChPT 
Wilson coefficients c^j and b k:i encode physics at the scale p 2 ~ A£ and the operators 
encode p 2 A£. 

At LO in the chiral expansion only one operator contributes in our analysis: 

Pf 



/np,a _ 
W k,0 ~ 



1 + (-l) k Tr [\ a LR n k tf - A^D fc £] , 



(23) 



where D k = (in ■ d) k , f p = f v M in the chiral limit. 0\f^ vanishes when k is odd due to charge 
conjugation and SU(3) symmetry. The normalization of 0\f^ is chosen such that c 0j o = 1- 
Note that the operator 

Pf 



Tr [Xl R ^ (D*E) St - \ R J1 (D fe St) E] 



is also LO but not independent of O^. Since □ fc (S'l'S) = 0, 

= (□ fe £ t )£ + £ t (D fc £) + ... , 



(24) 



(25) 



where the ellipse denotes (n h m £t)(D m E) terms that only contribute for matrix elements 
with more than one meson. Thus, Eq. (25) allows us to move factors of d k onto a neighboring 
£ and reduce <9£J to <9£J. 

At NLO there are two independent operators: 

Pf 



~rP' a — 
u k,i — 



O 



Pf 



k,2 



Tr 



1 + (-l) fc j Tr [m,Et + EmJ] Tr [A^tfE* - A^D fe £] , 
A^HED fc Et + (-l) fc (D fe Et) EmJ} 



-X a RL {m^Ws + (-l) fc (□*£) tfm q } . 



(26) 



(27) 



All other NLO operators have derivatives on more than one E, or can be reduced to 
and Or' using the equations of motion. For instance, consider 

Pf, 



O 



k,3 
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-Tr 



X a LR {^m q D k ^ + (-l) k (n fe £t) m ? Et} 



-AV{Smtn fc E + (-l) fe (□*£) mjs}]. 



(28) 



— ® — 

FIG. 1: NLO loop diagrams, where here (g) denotes an insertion of O^'q, and the dashed lines are 
meson fields. 

The sum and difference O k ' 2 ± 0^' 3 contain factors of (Y)m q ± mjE) and (Em£ ± m q Y$). 
We can trade O k \ — O k \ for operators with derivatives on more than one E by using the 
equation of motion for E, 

Et(^) 2 E = -(i^Et)(i9 M E)+S (S t m ff -mjE) , (29) 

together with the analogous equation for EL These operators with derivatives on more than 
one E do not generate one-meson matrix elements at tree level and can be omitted from 
our analysis. Thus only U^ 2 + O k '^ contributes and for simplicity we trade this for O k2 . 
We can also insert more E or Et fields into Sjg and get Tr [\ a LR T^m q T) (□ fc E) E+ + . . .]. 
But this operator can be reduced to O k2 using Eq. (25). Finally, we can consider operators 
where the power suppression is generated by derivatives rather than a factor of m q . Boost 
invariance requires that these operators still have k factors of n M , so they will involve D fc just 
like the operators we have been considering. To get power suppression with derivatives we 
can either use (c^E^c^E) which has derivatives on more than one E, or Et(<9 M ) 2 E which 
can be traded for operators with m q s using Eq. (29). Therefore, the operators with m g 's in 
Eqs. (26), (27) suffice. 

At NLO chiral logarithms can be obtained from loop diagrams involving the LO operators 
as shown in Fig. 1. For k = the operator O^ is the pseudoscalar current whose Fig. 1 
graphs give the one- loop corrections to For odd k the one-loop graphs vanish due to 



the 



-l) k 



factor originated from C. For any even k > the diagrams have a term 
where all derivatives act on the outgoing meson line, and this gives the same corrections 
as for f^. The first diagram could have additional contributions from derivatives acting 
inside the loop but it is straightforward to show that these diagrams vanish identically since 
n 2 = 0, and that the same holds true for LO operators with derivatives on more than one 
E [8]. Thus, we have shown that all possible non-analytic corrections are contained in f v M 
at NLO. This is true for every moment, and so we conclude that the leading order SU(3) 



violation of 4> p M (x) is analytic in m q . 

At NLO the Q' k2 and the wave function renormalization counterterms all contribute: 

(M b \c k , zH 2 O^ + ]T hM» = -if P M t ab (n ■ P) k (z% ■ (30) 
i=i 

For M = n and K, f p M can be related to by relating ip^\ a ip and the derivative of 
^7^75 A a, in the operator level. Acting Z}^ on both sides of 

(M^h^^O) = -i/ftJVe*' , (31) 

then by the equations of motion and Eq.(12) with k — 0, one obtains 

2 

Jir J* 2m > 

/£ = /l=^- • (32) 
m + m s 

These relations should be reproduced to all orders in ChPT. Indeed, direct computations of 
the leading chiral logarithmic corrections of f p Ml and m 2 M confirm these results. There 
is no relation between and however, because 

iD^r] \uYl 5 u + ~drf^ h d - 2s 7 ^ 75 s|0) = -2(r]\m (u lb u + d^d) - 2m s s l5 s\0) , (33) 

which is not proportional to (j]\u r )^ ) u + d^d — 2s7ss|0) away from the SU(3) limit. 
By comparing the k — and k ^ cases in Eq.(30), we obtain 

(z% = <^ + 2{(l + (-l) fc )TrK]^(&M-ViC fc ,o) 

+Tr [m q (\ b \ a + (-l) fc A a A 6 )] (6 fc)2 - b , 2 c k , (l + (-l) fe ) /2) }, (34) 

where (z k ) P = (z h ) P M in the chiral limit and (z°) p M = 1 is used. 
For k = 2m + 1 (odd moments), the {z k ) P M structure yields 

( z 2m+l)P = (z 2m+ X = 0, 

(z 2m+ % + = (z 2m+ % = (m s -m)b 2 m+i,2, (35) 
(z 2m+ %- = (z 2m+1 )^ = -(m s -m)b 2m+1 ^ 

For k = 2m (even moments), the (z k ) P M structure yields 



(z 2m )l = {z 2m )l + 2ma p 2m + (2m + m s )(3 p 



2mi 



(z 2m f K = (z 2m ) P + (fn + m s )a p 2m + (2m + m s )(3 p 2m , (36) 
( ^ = (^>-+^±i^o| m + (2m + m s )/3L, 
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where a\ m = (b 2m ,2 - &o,2C 2m ,o) and $L = 2 ( b 2m,i - b 0tl c 2m ,o). By isospin symmetry and 
charge conjugation, the even moments of different pion states (or kaon states) are equal. 
Eq.(36) implies a Gell-Mann-Okubo-like relation: 

(z 2m >£ + 3(^ m )P = 4^ 2m >^. (37) 

The moment relations imply the LCDF relations in Eqs.(4),(7). They also imply useful 
relations among the frequently used Gegenbauer moments, defined by 

2 



°" 2n + 



-jdxCl!\2x-l)<p p M {x) , (38) 



with a^f' p = 1. Here Cr/ 2 (z) denote the Gegenbauer polynomials which are even (odd) 
functions of z when n is even (odd). 



W=J2 a n' PC n 2 ( 2X - 1 )- (39) 



i=0 



Eqs.(35)-(37) imply that 



Km P = aZ + SaZ , (40) 

a 2m+l = a 2m+l = J (41) 

K°,P _ K+,P _ K*,P _ K-,P (An) 

a 2m+l ~ a 2m+l ~ a 2m+l ~ a 2m+l ■ K^J 



B. The <ff M Analysis 



Following the analogous procedures, the twist-3 operator (0£' a ) MJy can be decomposed 
into 

(orr = ra - (o a k r RL m 

with 

(OTYlr = i> L ^\ a LR \in-*D] k+1 i; R (44) 
and similarly for {0^ a ) R v L . Under charge conjugation, 

C -\0l' a y v C = {-lf{O a k ' a Y u (45) 

if we demand \ a LR ^ RL transform according to Eq. (19). At NLO in the chiral expansion, 0% 
is matched to LO and NLO hadronic operators O a k and O^: 

(orr - E c U°zr + E b u°Zr ■ m 



Again, at LO there is only one hadronic operator 



24 



1 + (-l) fe Tr [\ a LR (n [v & i] n k tf) - \ a RL (n [ "^n fc E)] , 



where n^d^ = (n^<9 M — n^d u ). At NLO, there are two independent operators, 



■a,a\ fiu 



rf 

24 



1 + (-l) fc j Tr [m q Yt + Sm[] Tr [\ a LR n [v d^U k T^ - \ a RL n [v d A U k 



(0Zr = ^Tr[A^{mtSn^n^t + (-l) fe (n^tfst)Smt} 
-\ a RL {m q ^n [u d^ ] n k J: + (-l) fe (n^tfE) E f m 9 } 

These yield the following relations. For /c = 2m + 1 (odd moments), 



{z 2m + l ): 



K+ 



(^ m+1 )^o = K-m)^ m+li2 , 



2to+1\<t 



For k = 2m (even moments), 



{z 2m+ %„ = -(m s -M)b° m+h2 . 



(z 2m )K = {z 2m )l + {m + m s )al m + {2m + m s )^ 
(z 2m ); = (z 2m )l + + 3 Ams) a L + (2m + 
where a£ m = (6^ m 2 - b^ mfi ) and = 2(^ m>1 - &o,i c 2m,o)- These im p!y and 

(z 2n x + 3(z 2m ); = 4(z 2m )^. 

In terms of the Gegenbauer polynomials 

oo 

<PUx) = 6x(l-x)Y,^C 3 J 2 (2x-l) , 
such that the Gegenbauer moments are 



i=0 



ai L " „ 4n + 6 „ ld.r(« 2 (2.r-\)<A,(.r) , 



n. 



6 + 9n+3n 2 



o 
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with a™' a = 1. Here Cl /2 {z) is an even (odd) function of z when n is even (odd), thus we 
have 



A K,<7 



u 2m ^ ou 2m 5 



l 2m+l ~~ "2m+l 



o, 



*2m+l 



A"+,cr 
a 2m+l _ 



~ a 2m+l 



—a 



K-,a 
2m+l 



(55) 
(56) 
(57) 



The normalization is related to for M = ir and K. By contracting 

M fe |^(x)a^ 75 A a i Z)H^(x)|o) = --/m^ 6 - f 9 m ) e ip ' x 

I 3 

with g va and making use of {M b ^(x)^X a il){x) |0) = f M 5 ab m 2 M p^e ip - x , one obtains 



(58) 



i-C-V 



m + m< 



m + m £ 
m K 



(59) 



Direct computations of the chiral logarithms of and /j£ also confirm these results. Again, 
f° and f£ are not related unless in the SU(3) limit. 



IV. THREE PARTON LCDFS 



The three-parton LCDF Tm for meson M is defined by [34]: 



= i 



{0\q'(yn) [yn,vyn]a fll/ >y 5 gG (Tp (vyn)\ a [vyn, -yn]q(-yn)\M b (p)) 
ifli5 ab [PaP[^9u] P - (o- <-> p)] T M (v, yp-n) , 



(60) 



where g is the strong coupling constant and G ap is the gluon field strength tensor, 
T M (v,yp-n) = J dae-* yp - n ( a '- a+va ^T M (a',a,a g ) , 

and 

Jda = J dada'da g 5 (1 — a' — a — a g ) , 



(61) 



(62) 
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and where a', a and a g are momentum fractions carried by q', q and the gluon. For simplicity, 
we use the gauge n • A — 0, such that the Wilson lines [yn,vyn] = [vyn, —yn] = 1. After 
the operator product expansion, the resulting twist-3 operator is 

°k°i = Yl ? { in ' d ) <VM (in - d ) A a (m • «9 ) g , (63) 



where we have suppressed the p, u, a, and p indexes on the left hand side. Note that O fc '° 
transforms in the same way as (0£' a ) A " y under chiral transformation and 

C-'O^C = (-l) fe Og . (64) 

Thus the matrix element 

(0|Oj a |M 6 ) = if M 5 ab [p aP]il g v]p -(a~ p)] (n ■ p) k+1 J da (a' - af oP g T{a', a, a g ) (65) 
implies 

%( v ) ot,OL g ) = T w ( v ) (a, a', a g ) , (66) 

T K + (a 1 , a, a g ) — T K + (a, a', a g ) = —T K - (a' , a, a g ) + T K - (a, a', a g ) oc (m s — m) , (67) 

2 [T K + (a 1 , a, a g ) + T K - (a 1 , a, a g )] = %(a', a, a g ) + 3T v (a', a, a g ) , (68) 

together with T K + = T K a and T K - = T-^o by isospin. Eq.(66) was also obtained in Refs. 
[33, 34]. 

The fact that Oj", (O u > a y u and O p Q ' a transform in the same way under chiral transfor- 
mation implies that the normalization constants /j^, and j v M receive the same leading 
logarithmic corrections. This is true between f v M and in Eqs. (32) and (59) due to the 
similar structures of the corresponding LO operators in Eqs. (23) and (47). Also, at NLO 
the ratio 

Rm = §- (69) 

JM 

satisfies 

Rk+ — Rk- ( m s — m) , (70) 

and 

AR K = R n + 3R V . (71) 
Analogous relations also exist for the ratio /m/ /m a ^ NLO. 
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As pointed out in Ref. [34] (see also [30, 33]), the three twist-3 LCDFs are related by 
the equations of motion. This implies in ChPT, the counterterms associated with the three 
LCDFs are also related. It can be shown straightforwardly that the SU(3) relations presented 
in this work are consistent with the constraints imposed by the equations of motion. 

V. LATTICE QCD CHIRAL EXTRAPOLATIONS 

The first few moments of the twist-2 and twist-3 LCDFs can be calculated using lattice 
QCD. Currently only the second twist-2 moment of pion (z 2 )^ is computed [24-26]. The 
latest results on (z 2 )^ have converged to (z 2 ) MS (// = 2.67 GeV 2 ) = 0.280(49)lffi [25] 
and (z 2 ) MS (fi 2 = 5 GeV 2 ) = 0.281(28) [26] with the lightest m n to be 490 MeV and 550 
MeV, respectively. In these two calculations, and many others, chiral extrapolations to the 
physical point = 140 MeV were performed. Thus the model-independent results on the 
m q dependence of LCDFs obtained in Ref. [7] and this work are valuable to reduce the 
systematic errors in the extrapolations. 

A very interesting result we found both in the twist-2 and twist-3 LCDFs is that, in the 
continuum, all the moments of those distributions are analytic in m q at 0(m q ). Remarkably, 
this nice feature is retained in a lattice theory with (partially) quenched, finite volume and 
the leading finite lattice spacing effects included. It is easy to see why. In all the cases, 
the leading logarithms are generated by one loop diagrams with insertion of the LO C^o 
operator [with the trace replaced by super-trace when (partially) quenched]. These diagrams 
correct Cfc j0 by a universal (k independent) logarithm which can then be absorbed into the 
normalization constant. Thus the moments are analytic in m q . This feature, however, is 
only true at 0(m q ). At 0(m 2 ), the m 2 In 2 m q corrections to c^o are generated by two loop 
diagrams with insertion of the LO Ok t o operator. These corrections are k independent and 
can be absorbed into the normalization factors. The m 2 In m q corrections, however, can arise 
from one-loop diagrams with insertion of the NLO Ok,i and O fc 2 operators which depend on 
bk t i and b k 2 and correct c^o in a A; dependent way, thus the m 2 \nm q corrections cannot be 
absorbed by the normalization constants. Generalizing the above argument to 0(m % ), the 
m % q In* m q corrections can be absorbed by the normalization constants, thus asymptotically, 
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the LCDFs has the quark mass dependence 1 

oo i— 1 
i=l j=0 

where a 00 and Qj{j £1X6 independent of m q . 

For different types of lattice simulations, the leading chiral corrections from the one 
loop diagrams can still be absorbed into the normalization constants, such that the LCDFs 
remain analytic at NLO. Different types of simulations give different normalization constants 
and different analytic contributions to the LCDFs. For example, the (2m + m s ) term in 
Eq.(36) is originated from Tr[m g ] in Eq.(34). In the (partially) quenched theory, Tr[m g ] is 
replaced by Str[m g ] which is zero in the quenched theory [35] and is {Tm sea + m SySea ) in the 
partially quenched [36] theory, with m qySea denoting sea quark masses. It is easy to check 
that Eq.(37) still holds in spite of (partially) quenching. As for finite volume effects, the 
counterterm structures are the same as those in the infinite volume calculations since the 
effects are infrared in origin. Thus in simulations with twisted [37-39] or partially twisted 
[40, 41] boundary conditions for fermions, the leading finite volume effects are absorbed by 
the normalization constant and the moments are all identical to those in the infinite volume 
calculations. For lattice fermions whose chiral symmetry is broken by finite lattice spacing 
a, new counterterms in powers of a need to be added to account for the leading finite a 
effects. Since these counterterms are SU(3) symmetric, Gell-Mann-Okubo-like relations like 
Eq.(37) will still be true. The same is true for mixed action simulations with chiral valance 
quarks but Wilson or staggered sea quarks. 

1 Consider a generic diagram giving an 0(e d ) contribution to the twist-2 matrix element {M (p)\Of? \0) , 
where e ~ p ~ m^. If the diagram has a £)(e fc + 1 + 2n ) pionic operator insertion, / internal propagators, L 
loops, and Vi vertex of order e di , then d = fc + l + 2n — 2/ + 4L + ^\ V^d,. After using the topological 
identity L = I — J2i V an d removing the e k+1 kinematical factor, this diagram gives a 0(e 2h ) contribution 
to the moment (z k ) M with h = n + L + J2i Vi (di/2 — 1). Since n > and d% > 2, h > L. In general, 
the diagram has the contribution ^2f =0 biirig \n m q , with the power of logarithm not bigger than L and 
m q oc ^/jtitt = 0(\fe). The In' 1 m q term can only come from the diagram with an insertion of the 
LO pionic operator (n = 0) and can be absorbed into the normalization. The arguments can be easily 
generalized to the twist-3 cases. 
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VI. COMPARISON WITH LIGHT CONE SUM RULE RESULTS 



The first few Gegenbauer moments for the octet mesons have been computed using light 
cone sum rules [27-32]. The twist-2 moments a% I ' P (defined in the same way as a^ ,a in 
Eq.(54)) were determined in different references. Quite different values for al' P (/j, = 1 
GeV), 0.44 [27] and 0.20-0.25 [28], are obtained using light cone sum rules. After combining 
these values with various constraints from experimental data [42, 43], a reasonable range is 
assigned [43]: 

< a 2 ' P (l GeV) < 0.3 . (73) 
For moments of K and 7], light cone sum rules of Refs. [32] and [30] yield 

af' P (l GeV) = 0.16 , a v / '(1 GeV) = 0.2 , (74) 

respectively. Since the range of a 2 P is big, the results in Eqs.(73),(74) are still consistent 
with the ChPT relation found in Ref. [7]; 



3a 



2 i , /?V™2 



1 + 0(m ) . (75) 



However, this equation could provide a tight constraint among a% r,P since one expects the 
0(m 2 ) correction to be 10% at most. 

As for af ~' P , recent results yield af ~' P = 0.05 ± 0.02 [44], 0.10 ± 0.12 [45] and 0.050 ± 
0.025 [46]. The positive sign corresponds to the s-quark in a K~ meson carrying a bigger 
momentum fraction than the w-quark. 

In Ref. [32], twist-3 Gegenbauer moments are computed (see also [47]): 

M,p o M M 27 2 81 2 M,P /ry r \ 

a 4 •* = -3r7 3 uj 3 - —p M - ^Pm° 2 , ( 76 ) 

M,a r M ^ M M 7 2 ^2 M,P 

a 2 = 5% - -% ^3 - t^Pm - -^PmO>2 ■ 

The numerical values of r/^ 1 and uj^ 1 are independent of M, while p 2 M oc m 2 M = 0{m q ) and 
a^' P = a P (l + 0(m q )) with a P = a^' P in the chiral limit. Eq.(76) satisfies the relations 
Eqs.(40),(55) and, hence, is consistent with chiral symmetry. 



15 



VII. CONCLUSIONS 



Using chiral symmetry we have investigated the leading SU(3) violation in the complete 
set of twist-3 light-cone distribution functions of the pion, kaon, and eta, including the 
two-parton distributions (f> p M , (ff M and the three-parton distribution <$? M . It has been shown 
that terms non-analytic in the quark masses do not affect the shape, and only appear in 
the normalization constants. Predictive power is retained including the leading analytic 
m q operators. With the symmetry violating corrections we derive useful model-independent 
relations between (f) P r ' a , <pP' a , 4> P R+ R0) and 4> p ^ K _. We have also commented on the calculation 
of the moments of these distributions using lattice method and light-cone sum rules. 

We thank Will Detmold, Dan Pirjol and Iain Stewart for useful comments on the 
manuscript. Support of NSC and NCTS of ROC is gratefully acknowledged. 
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